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Besidea eatablishlng the uae of an Interactive systea to obtain 
eatloates of aoll oolature* the preaent reaearch denonatratea that there la 

aoae Juatlflcatlon for attempting to determine pattema of molature 

\ 

availability on the regional scale. The function of a boundary layer model 
la to derive the molature availability parameter from the aeaaured ground 
temperature reaponae. Although there la some conceptual Indetermlnancy 
In the model concerning the real nature of the surface, the moisture 
availability la probably a purer representation of the surface moisture, and 
certally more Intrinsic to the ground surface, than la the temperature 
response Itself. 

More than conceptual problems stand In the way of deriving meaningful 
aoll moisture values. It Is doubtful If any satisfactory representation of 
the vegetation canopy can be made for general purposes, although It Is 
conceivable that small-scale remote measurements of the blackbody tempera- 
ture may lend themselves favorably to a more precise parameterization of the 
surface canopy where the canopy type and height Is known. More field mea- 
surements are necessary to Improve the vegetation parameterization. Another 
serious Impediment to satellite measurement techniques Is the Interference 
by cloud and moisture. Cirrus or scattered cumulus can go undetected 
by conventional analysis but the presence of such obstructions can lead to 
anomalously cool or Irregular temperature patterns. Moisture Inhomo- 
geneltles can also produce similar, though less drastic, anomalies. Because 
of clouds, higher resolution satellles, such as HCMM, afford less problems 
In data Interpretation. Indeed, our results seem to de^ownstrate some 
superiority of HCMM over GOES In that respect. 
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The feet thet cloud le really aore prevalent than one would Initially 
suppose froa casual Inapectlon of satellite photos, leads to a problea In 
continuity of aeaaureaent. It Is not easy Co find two or three clear-sky 
satellite aeaaureaents over Che saae area on a given day, a necessary 
arrangeaent for deteralnlng a soil aolsture estlaate. Since the purpose of 
this research originally was to derive patterns of aolsture availability 
which are needed In atmospheric models to obtain realistic patterns of 
surface sensible heating. It may be that an archive of moisture availability 
patterns can be generated from satellite data for a particular region. 

For agricultural use the precise meaning of the moisture availability 
may be slightly ambiguous except as a rough index of aolsture stress, but 
for numerical atmospheric models the Importance of the moisture availability 
is that it will yield something approximating the correct surface heating 
patterns, which are so necessary in predlcitng the evolution of weather 
systems. Finally, we suggest that existing HCMM satellite data constitutes 
a great wealth of information for further model testing and method 
improvement. 
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Abitract 

An interacclve method for obtaining patterns of moisture availability 
(and net evaporation) from satellite infrared measurements has been 
developed at The Pennsylvania State University. The method employs the 
boundary layer model of Carlson (Carlson et el.. 1981) » and a variety of 
image processing routines executed with the aid of a minicomputer. 

To test the method vith regard to regional scale moisture analyses, two 
case studies were chosen for the availability of HCMM data and because of 
the presence of a large horizontal gradient in antecedent precipitation and 
crop moisture index. Results showed some correlation in both cases between 
antecedent precipitation and derived moisture availability. Apparently, 
regional-scale moisture avalability patterns can be determined with some 
degree of fidelity but the values themselves may be useful only In the 
relative sense and be significant to within plus or minus one category of 
dryness over a range of 4 or 5 categories between absolutely dry and field 
saturation. Interpretation of surface moisture availability patterns 
undoubtedly must vary according to the nature of the surface (e.g. 
vegetation type) but our preliminary results suggest that the derived 
moisture values correlate best with longerrterm precipitation totals, 
suggesting that the Infrared temperatures respond more sensitively to a 
relatively deep substrate layer. 
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I . Incroductiofi 

Drawing on ch« prallBinary but soaewhac suggestlva results of Kocin 
(1979) and Cooper (1981)* which indicate that soae Masure of soil ooisture 
can be determined remotely over vegetated canopies using thermal Infrared 
measurements » we attempted to use HCNM infrared measurements to determine 
regional*^cale patterns of moisture availability in two case studies. These 
case studies were chosen on the basis of available HCMM imagery (a day-night 
image pair within a 36h interval), absence of clouds and the presence of a 
large horisontal variation of precipitation and crop moisture index (CMI) 
across a region. We reasoned that if moisture availability patterns derived 
from the thermal method described by Carlson et al. (1981; henceforth 
referred to CD) are to prove meaningful then they should, in some way, 
reflect the distribution of antecedent precipitation and that the 
relationship between moisture availability and precipitation should be most 
obvious over regions containing sharp variations in aridity, as measured by 
the crop moisture index. Accordingly, we chose one case study for August 
1978 (August 22-23) over southern Indiana, Illinois and northern Kentucky 
and another for July 26-27, 1978 over eastern Kansas. In both of these 
cases the crop moisture index varied from about - 2 to -f 1 across the 
domain. 

The main point of departure between this investigation and those 
reported by Carlson et el. (1981) or Carlson and DiCrlstofaro (1981) for 
urban areas and by Kocin (1979) is the scale. We were specifically inter- 
ested in the regional scale (a few hundreds of kilometers on a side, 
sometimes referred to as mesoscale) because we wanted to learn whether the 
CD could be applied to determine useful patterns of moisture availability 
and thermal inertia on a scale commensurate with existing mesoscale 
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precipitation natworka, ao that «e would be able to make use of the antece- 
dent precipitation as a kind of ground truth. In addition, we felt that the 
application of ■atelllte-derlved paranetera to ataoapherlc nunerlcal 
modeling and prediction could more easily be demonstrated If these patterns 
were constructed over mssoscale regions because of current Interest In 
weather prediction on that scale. A third reason for expanding the scale 
from that In earlier studies was to compare patterns derived from HCMM 
measurements, with those from the geostationary satellite (GOES), whose 
resolution is about 4 km. 

Besides, the case study approach an Important objective In this 
research has been to further develop an Interactive approach for determining 
soil moisture patterns with the idea that we will achieve a nearly 
operational capability using a olnlcomputer-lmage processor system available 
In the Department of Meteorology. 

The essence of our work during the past year Is summarised by Polansky 
(1982) in his M.S. thesis, a copy of which Is attached as Appendix II for 
convenience of reference. All figures, unless otherwise noted, will refer 
to those In Polansky 's thesis. 

2 . Results 

2.1 Case studies: Kansas and Indiana 

Polansky reworked two case studies Initiated by the principal Investi- 
gator. His research was divided Into four general areas: model sensi- 

tivity, error analysis, analysis of moisture, and system development; the 
last was omitted from the thesis. This report %rlll confine Itself to a 
brief account of the results and error analyses discussed by Polansky. 

In order to derive patterns from HCMM (resolution 0.5 ka) appropriate 
to the scale of GOES, HCMM pexels were combined to form a larger pixel about 
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2 km on a side consisting of about 20 original pixels; the 128x128 pixel ^ 

working area thus contained an area about 250 km on a side. 

In the Indiana caae, a largely rural, vegetated area (Fig. 4.1), was 
experiencing a large horlsontai variation in crop moisture Index and 
antecedent (3 week total) precipitation with serious drought conditions (crop 
moisture contour value of 1 omitted) developing over southern Illinois 
where less that 0.25 Inches of precipitation had fallen in 3 weeks. Despite 
abundant cloud cover, which largely ruined the Impact of the results. It Is 
apparent in Fig. 4.3A that the driest region was also the wariMst. In the 
moisture availability and total dally evaporation patterns derived from CD 
(Figs. 4.4, 4.6), the lowest values of N «0.25) appear over the driest 
region. Thermal Inertia (P) values also tended to be a little lower over 
the ralnfall-deflclt region (Fig. 4.5). 

The Kansas case proved to be better for Illustrating a relationship 
between crop moisture Index and satellite-derived moisture patterns because 
there was less cloud and a greater horizontal variation In crop moisture 
Index (Fig. 4.13B, 4.14) over a region consisting largely of unlrrigated 
crop and grazing land (Fig. 4.12B). Warmest temperatures during the day 
(Fig. 4.15A) coincided closely with the region of low crop moisture Index 
and light precipitation amounts. However, the temperatures appear somewhat 
too high, at least as compared with those of GOES (Fig. 4.21B). Moisture 
availability (Fig. 4.16) and total evaporation (Fig. 4.19) were also 
relatively low and the heat flux relatively high (Fig. 4.18) over the arid 
southwestern portion of the Image area. GOES moisture heat flux malyses 
were rather similar to those of HCMM, though differing In detail. 

It Is evident that In both cases low moisture availability and high 
surface sensible heat flux tended to coincide with regions where 
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prtcipltatlon aaounts w«ro leas than 0.25 to 0.50 Inches (and a crop 
aolsture index of less than - 1) and that relatively high iM>lsture avail- 
ability and lov surface sensible heat fluxes tended to be located where the 
three-week precipitation totals were in excess of about 1.00 Inches (crop 
nolsture Index of greater than zero). The correspondence between the 
aridity Indices (antecedent precipitation or crop moisture) Is obviously 
very rough. A more systematic view of such a relationship Is presented In 
Pig. 1 of this report In which the antecedent precipitation, as defined by 
Blanchard et al. (1981; unpublished manuscript), Is plotted against moisture 
availability. In the figure the data was read from corresponding grid 
points on both the rainfall and moisture availability charts. Despite con- 
siderable scatter the graphs exhibit some systematic behavior, as would 
be expected from a qualitative examination of the figures. Interestingly, 
a better relationship was obtained for the antecedent precipitation 
weighting factor (k) of 0.93 than of 0.78, suggesting that the moisture 
availability (and hence the surface response measured by satellite) Is more 
sensitive to long-term precipitation trends, and hence, to a relatively deep 
rooting depth. For comparison, a similar analysis Is presented for Kocln's 
Missouri Watershed analyses (Fig. 2) except that the precipitation In the 
figure simply refers to a cumulative three-week total. It Is to be 
expected, however, that the relationship between moisture availability and 
antecedent precipitation will vary with surface type. There Is apparently 
some relationship between precipitation and moisture availability, at least 
In the June case (Fig. 2a), but the Septemter analysis (Fig. 2b) shows very 
little correlation between the derived moisture availability and antecedent 
precipitation. For comparison. Cooper's results of nolsture availability 
versus antecedent precipitation (weighting factor 0.93) Is presented In 
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Fig. 1; tha bar langth danotaa uncarcalnty dua to tha ranga of aurfaca f 

camparacuraa oaaaurad ac dlffarlng anglaa ovar diffarlng crop typaa uaing 
a handmaid radioratar. 

Tha ralacivaly waak corralatlona avldant la Flga. 1-2^ la not aacaa- 
aarlly a cauaa for dlaeoaragaoMnC. Antacadant praclpltatlon la not a trua 
aaaaura of aoll oMlatura baeauaa tha lattar can ba aa aaally Influancad by 
tha praaanca of undarwatar acqulflara, varlatlona In watar tabla laval, 
dlffarancaa In vagatatlon typa and aga» runoff, and artificial watering 
(a.g. Irrigation) aa by rainfall Itaalf. Although tha crop oolatura Idax 
rapraaanta an attanpt to Incorporate runoff and evaporation In dataralnlng 
an affective maaaura of aurfaca oolatura atreaa. It la datarmlnad for 
avaragea ovar relatively large araaa. Ware tha boundary layer model and 
aurfaca temparatura meaauraoanta perfect, one would naverthaleaa anticipate 
a high degree of acattar between rainfall and nolature availability, even on 
tha acale of a HCMM pixel. 

2.2 Errora and lloltaciona 

The reaulta of our two caaa atudlea, aunoarlzed In Flga. 1>2, are 
alnlnar to that of Harlan (1981) who plotted antecedent praclpltatlon veraua 
day-night teoperatura dlffarancea and found a ayatenatlc Inverae relatlon- 
ahlp, albeit one affected by a conalderable acatter of data. There la an 
unmlatakeable relationahlp between effective aoll molature, vaguely 
repreaented by antecedent rainfall, and tha derived oolatura avallabllty, 
although It la not poaalbla to define praclaaly tha aurfaca or aubatrata 
(rooting) depth appropriate to tha Inferred value of oolatura availability. 
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Tha beat correlation between oolatura availability and antacadant 
precipitation (r 2 - 0.4) waa for tha HCMM Kanaaa data. 


Quit* probably tha lataauramone la Intiaataly llnkad In aach pixal to tha 
typa, haight and aga of tha vagatatlon canopy and a boat of othar aubatrata 
factora uhlch govam tha ratantion of watar. An ataoapharie aodal, auch aa 
tha <ma uaad to obtain Polanaky'a raaulta, aarvaa only to raaova that 
variability In tha aatalllta aaaauraMnta Which la dua to fluctuatlona In 
tha ataoaphara* Although varloua ataoapharlc laodala othar than tha CM (a.g. 
•odala of Prlca, 1982; Wataon and Miller, 1981; Pratt and Ellyatt, 1979; 
Roaaaa, 1978), would likely yield coaparabla ceaulta, there la probably an 
Irraducaabla lower Halt of error of t 1*>2*C, below which no aodel la 
capable of attaining. Thera are many raaaona for thla Inherent error auch 
aa lack of underatandlng and formulation of tha correct boundary layer 
phyalca, uncertainty about the nature of aurface and aubatrate, 
non-advectlve affecta In the atooaphere, errora In Input data (l.e. In 
teaperature) due to clouda, Inatruaentatlon, and water vapor attenuation. 
There la al&o the problem In averaging a aeml -Infinite mimber of aurface 
Infrared emmltora to produce a alngle pixel maaaurement. Thua, errora may 
be dependent on both model and radiometer. Aalde from definable errora 
there la alao a nagging conceptual problem In Identifying that particular 
aurface to which a particular molature availability pertalna. 

Becauae of theae errora, there are only certain comblnatlona of 
a'itelllte Image tlmea aultable for uae. We have Inveatlgated three typea of 
image comblnatlona for determining aoll molature: a day-night Image pair, a 
aornlng-afte moon-night triplet and a momlng-afternoon pair. Only one of 
theae comblnatlona, tha flrat one, la appropriate to HCMM orblta but, aa It 
tuma out the HCMM overpaaa tlmea were virtually optimal for day-night Image 
palra. Nevertheleaa, even with the auperlor reaolutlon of HCMM It la likely 
that the derived molature availability la i«jbject to an error of i 0.1-0. 2 
and, conaequantly, molature availability pattaraa Juatlflably might be 
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rapr«s«at<d wiChouC fux lota of datall by four or flva quailcatlvo caco* 
gorlea: a.g. vary dry < 0.2); dry (M ■ 0.2-0. 4); narglnally wat (M • 

0.4-0. 6); wat (M ■ 0.6-0. 8) and vary vat (N ■ 0. 8-1.0). 

3. Conclualona 

Ragional acala aaaauraaanta of aoil nolatuta ualng aatallita Infrarad 
taaparatura aaaaurananta nay ba poaalbla but Cha raal validity of tha 
Boiatura availability valuaa can not ba aaaaaaad without ground truth naa- 
aurananta; So far auch diract aaaauranai?ta ara iavractical to obtain and 
tha only indiract avidanca of a ralationahip batwaan darivad noiatura 
availability and aoil aoistura ia inparfactly raflactad in cha pactarn of 
antacadant rainfall. A high-degraa of aci.cter in cheaa raaulca may ba a 
conaequanca ;)f aoil noiatura fluctuationa raaulcing fron factora othar than 
precipiCation» auch aa rooting dapth of the vagatation. runoff irrigation, 
ate. 

In two caaa aCudiaa low crop noiatura index and anall valuaa of 
antacadant precipitation correlated with low valuaa of noiatura 
availability. Drought conditiona vara clearly related to higher afternoon 
aurfaca tanparaturaa , which ara equivalent to lower valuaa of noiatura 
availability according to the boundary layer nodal. Tha noat aarioua 
problan, however, ia not tha nodal icaalf but tha conceptual franevork 
which ia linitad in tha face of an anornoualy conplax and ground or vege- 
tated aurfaca. Thua, while any raaaonabla atnoapharic nodal appaara to 
ba capable of ranoving external atnoapharic influancaa on tha aurfaca 
tanparatura variation and thereby provide a purer naaaura of aoil noiatura 
chan doaa cha bulk day-night tanparatura difference, there ranaina cha 
quaation:-- what actually doaa tha noiatura availability repraaanC? It nay 
ba, however, that ona nuec aattla in tha and for a noiatura availability 
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p«raMt«r which la abla to ropraaant tho raal toil ■olatura only Indiroetly 
and trhlch Itaalf say ba accurate to within plus or ainua one category in a 
range of 4 or S categoriee of soil dryness. 

Finally, in coapering GOBS with HCMf, patterns of aoisture availabllty 
seea to be aore seniitive to surface precipitation with HCMM, auggeating 
that the finer resolution snd generally superior quality of HCWt data is 
better suited to deriving soil aoisture estiaates. In any cane, we now 
possess the capability of producing aoisture svailsblilty asps over a saall 
scale or aesoscals region interactively in a aock operational node using the 
Penn State ainicoaputer syscea and iaage processor. 

A coaplete description of the Infrared asthod will be contained In a 
forthcoming paper by Carlson (198r). 


Appendix I 


List of HCHM scenes for which analyses were perfonsed on 
the coaputer coapatible tape data. 


Alls 185801,2 

Indiana 

day - 

August 23, 

1978 

All? 074803 

Indiana 

night - 

August 22, 

1978 


A0093 193601,2 

Kansas 

day - 

July 28, 

1978 

A0092 082303 

Kansas 

night - 

July 27, 

1978 
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Figure 


Figure 


Figure Captions 


1. Moisture availability versus antecedent precipitation (API In 
Inches; weighting factor k * 0*92) for two case studies using 
HCMM data for Kansas (July 27-28 » 1978; dots) and Indiana August 
22-23, 1978; crosses). Values represent corresponding 
measurements at resampled pixel points. T^e thin dotted line 
modifying the linear correlation, suggests an alternate 
relationship between M and API. 

2a. Same as Fig. 1 but for cumulative three-week precipitation 
amounts versus moisture availability for Kocln's (1979) 

June 9-10, 1978 HCMM Missouri Watershed analyses 


Figure 2b 


Same as Fig. 2a except for September 29, 1978 


M «rs. APr For HCMM KS/I//D 

i? iTfS 

, KS cxse. 


ORIGINAL page » 
OF POOR QUALfTf 


X rvo Mt 


/ f\PT- 0«XO **• 

: 0.31 


•/' 


- ^ ' 


CIO 

3 

C ^ tc Zi 

C 3 C 

•H < ^ 

•-» Cfi ^ 

a, 3 - 

< C -3 

w CO li 

0^ ^ tC 4J 

C •*- T3 c -u 

O 73 C ^ C 

^ 3 ^ 

w C 

JO Cfl -3 O C 

4J C 3. 

•H tt> 3 'J5 £ 

a. V) 1 ) 4j 

0 U CO W 3 

3 w 0 ^ 

L- O O U H 

a 3 73 

4J 4J 

4J C • 

c u CO 3 ce 

V O 3) w 

XI ^ O' a# e 

01 ^ W -p4 

0 ^ e- 0 

01 • 01 & 

4J O' 00 bt 

C • <M ^ 

3 0 1 CO 01 

01 K 

CO R CM 3 *<Ni 

3 *H D. 

QO JE >s 3 

W <-4 > X3 

3 Lri 3 3 

> O ^ ^ 

iJ v.y^ • 3, 

3s y 8 

AJ 3 3 3 3 

•«i4 <4^ 3 y 3 

^ CO 3 3 

«p<4 00 C 3 bi 

^330 
3 *H ^ U 4J 


•/ ^ 
I • 


U •• 3 
3 3-0 
U 3 
3 ^ X 

O S 
o e o 
r ^ X 


r 


Mv5. apt Missouri 

ORIGINAL PAGE IS 
OF POOR QUALmr 



Sa r¥» as Fig. 1 but for cumulative three-week precipitation 
amounts versus moisture availability for Kocln's (1979) 
June 9-10, 1978 HCMM Missouri Watershed analyses 
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Figure 2b, Same as Fig, 2a except for September 29 



